Dynamical approach to spectator fragmentation in 
Au+Au reactions at 35 MeV/A 
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The characteristics of fragment emission in peripheral 197 Au+ 19 ' Au collisions 35 MeV/A are stud- 
ied using the two clusterization approaches within framework of quantum molecular dynamics model. 
Our model calculations using minimum spanning tree (MST) algorithm and advanced clusterization 
method namely simulated annealing clusterization algorithm (SACA) showed that fragment struc- 
ture can be realized at an earlier time when spectators contribute significantly toward the fragment 
production even at such a low incident energy. Comparison of model predictions with experimental 
data reveals that SACA method can nicely reproduce the fragment charge yields and mean charge 
of the heaviest fragment. This reflects suitability of SACA method over conventional clusterization 
techniques to investigate spectator matter fragmentation in low energy domain. 

PACS numbers: 25.70.-z 25.70.Mn 24.10.Lx 



I. INTRODUCTION 

The study of heavy-ion (HI) reactions at intermedi- 
ate energies provides an important platform to probe the 
highly non-equilibrium environment produced in the re- 
action zone as well as properties of excited fragments 
evolved from the spectator zone [H-Q- The spectator 
matter fragmentation at relativistic bombarding energies 
is also characterized by a rise and fall pattern [J-H, 
These experiments performed on ALADiN set-up mainly 
focussed on the liquid-gas phase transitions and uni- 
versality behavior observed in the fragment-emission at 
incident energies > 400 MeV/A [2|, |5j. The clusteri- 
zation approach based upon minimization of fragments' 
total energy dubbed as simulated annealing clusterization 
algorithm (SACA) is reported to explain this rise and fall 
trend in the multiplicity of IMFs with collision geometry 
quite accurately at relativistic bombarding energies \J\. 
Conventional clusterization approach based upon spatial 
correlation among nucleons, however, failed completely 
to explain higher IMF yields at large impact parameters 
01- 

The reaction dynamics and associated non-equilibrium 
aspects in low-energy HI collisions are, however, still 
poorly understood phenomena 0, Hoj |. The important 
question associated with this domain is whether nuclear 
system can reach thermal equilibrium before break-up 
or not [11]. Most of nucleon-nucleon collisions are Pauli 
blocked and mean fiel d g overns the reaction dynamics 
in low energy domain [r3 - [l4l ]. In case of peripheral col- 
lisions, we have projectile-like and target-like remnants 
and fusion-fission events dominate the scenario. The fu- 
sion events disappear as the incident energy increases 
marked by the onset of multifragmentation. In the low 
energy range {i.e. between 40 and 70 MeV/A), a signif- 
icant fraction of the intermediate mass fragments orig- 
inates from mid- velocity region (s}. This type of pre- 
equilibrium emission has been conjectured as 'extended 
neck' emission due to the dynamical fluctuations that in- 
crease with the incident energy. Many experiments have 



indicated that binary dissipative collisions (BDC) dom- 
inate the scenario in such low energy HI collisions (l5j . 
As far as decays from quasiprojectile and quasitarget are 
concerned, molecular dynamics approaches coupled with 
conventional clusterization algorithm mayn't reproduce 
fragment charge yields accurately [l|| . 

In the present paper, we aim to see whether SACA 
method can describe spectator matter fragmentation in 
low energy domain or not. The confrontation of the the- 
oretical predictions employing advanced clusterization 
technique with experimental data can be of importance 
to understand the physical scenario behind cluster pro- 
duction mechanism in low-energy HI collisions. We shall 
compute the fragment observables for peripheral Au(35 
MeV/A) + Au collisions employing MST and SACA clus- 
terization subroutines. A comparison with experimen- 
tal results recently obtained by Multics-Miniball group 
[Til [HI is also attempted to explore the applicability of 
SACA method. Section [TT] describes the main features of 
the QMD model along with simulated annealing cluster- 
ization algorithm (SACA). Our results are discussed in 
section HTT1 and summarized in section HVl 



II. 



THE MODEL 



The quantum molecular dynamics model is yl-body 
transport theory that incorporates the quantum features 
of Pauli blocking and stochastic nucleon-nucleon (n-n) 
scattering [TJ, [l9| • Each nucleon in the colliding system 
is represented by a gaussian wave packet as 
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Mean position Ti(t) and mean momentum Pi(t) are the 
two time dependent parameters. The gaussian width 
has fixed value v / L=1.08 fm. The centers of these 
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Gaussian wave packets propagate in coordinate (7^.3) 
and momentum (V3) space according to the classical 
equations of motion: 



Pi = 



9(H) . 
dr.; 



r, = 



d(H) 
dpi 



(2) 



The Hamiltonian T-L appearing in Eq. [2] has contribu- 
tion from the local Skyrme-type, Yukawa and effective 
Coulomb interactions [14j • Since QMD model gives just 
the phase space of nucleons, one needs secondary al- 
gorithm to clusterize the phase space. The standard 
clusterization approach namely minimum spanning tree 
(MST) procedure assumes that the correlating nucleons 
belong to the same fragment if their inter-nucleon dis- 
tance I Yi — Yj I is smaller than 4 fm. This approach is 
successful when the system is dilute and clusters are well 
separated in IZ3 space. 



A. The SACA formalism 

This sophisticated clusterization approach allows early 
identification of fragments before these are well separated 
in coordinate space. The SACA method works on the 
principle of energy minimization of fragmenting system. 
The pre-clusters obtained with the MST method are sub- 
jected to a binding energy check [20l[2lj: 
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< —Ebind, (3) 



with Ebind = 4.0 MeV if Nf > 3 and E^nd = otherwise. 
In Eq. ([3]) , Nf is the number of nucleons in a fragment 
and P^" 1 ' is the center-of-mass momentum of the frag- 
ment. The requirement of a minimum binding energy 
excludes the loosely bound fragments which will decay 
at later stage. To look for the most bound configuration 
(MBC), we start from a random configuration which is 
chosen by dividing whole system into few fragments. The 
energy of each cluster is calculated by summing over all 
the nucleons present in that cluster using Eq. ([3]). 

Let the total energy of a configuration k be Ek(= 
Ei^YfCi)) where Nf is the number of nucleons in a frag- 
ment and Q is the energy per nucleon of that fragment. 
Suppose a new configuration k (which is obtained by (a) 
transferring a nucleon from randomly chosen fragment to 
another fragment or by (b) setting a nucleon free, or by 
(c) absorbing a free nucleon into a fragment) has a to- 
tal energy Eg. If the difference between the old and 
new configuration AE(= E k > — Ek) is negative, the new 
configuration is always accepted. If not, the new configu- 
ration k may nevertheless be accepted with a probability 
of exp(—AE/v), where v is called the control parameter. 



This procedure is known as Metropolis algorithm. The 
control parameter is decreased in small steps. This algo- 
rithm will yield eventually the most bound configuration 
(MBC). 

Since this combination of a Metropolis algorithm with 
slowly decreasing control parameter v is known as sim- 
ulated annealing, so our approach is dubbed as sim- 
ulated annealing clusterization algorithm (SACA). For 
the present calculations we have employed improvised 
version of SACA method in which fragments are con- 
fronted against realistic binding energies instead of con- 
stant binding energy check of -4.0 MeV/nucleon (2^. It 
may be mentioned that this modification has nearly no 
effect on final fragment yields and it explains the AL- 
ADiN multifragmentation data @ at relativistic energies 
quite nicely. 



III. CALCULATIONS AND COMPARISON 
WITH DATA 

We simulate the reactions of Au (35 MeV/A)+Au 
at six peripheral geometries using the MST and SACA 
methods. A soft equation of state (k — 200 MeV) along 
with standard energy dependent n-n cross section was 
used for the simulation of HI reactions. 

Figure 1 displays the average density reached in the re- 
action along with the evolution of the heaviest fragment 
A max and multiplicity of clusters with mass A > 5 fol- 
lowed till 300 fm/c. One can see that as the nucleon den- 
sity saturates after violent phase, SACA method is able 
to identify asymptotic size of the heaviest fragment A max 
around 100 fm/c. The QMD+MST approach, on other 
hand, fails to detect the heaviest fragment A max even at 
300 fm/c. The heavier A max detected with MST method 
at 100 fm/c implies that it still assumes smaller fragments 
as constituents of A max , being very closely spaced. We 
shall show later on in 3-D space that this A max in actual 
consists of bunch of smaller clusters. The higher multi- 
plicity of fragments with mass A>5 obtained using SACA 
method (see figure lc) clearly shows that A max detected 
in the MST approach is actually a bunch of smaller frag- 
ments. As shown in Ref. [21| . different clusterization 
approaches converge to same configuration at 1000 fm/c. 
It is, however improper to follow the reaction till 1000 
fm/c with semi-classical theory, when nuclei are found to 
emit fraction of nucleons even after 100 fm/c. 

In figure [U we explore the sensitivity of fragment 
charge distribution towards clusterization algorithms. 
We display here the calculations using MST and SACA 
approaches for unfiltered events along with experimen- 
tal charge yields (shown as solid circles) obtained from 
the decay of the quasiprojectile (QP) [17]. Top and bot- 
tom panels show the comparison of model predictions 
with experimental data in the impact parameter inter- 
vals 0.5 < b/bm ax < 0.6 and 0.6 < b/b max < 0.7, respec- 
tively. The MST method clearly predicts larger produc- 
tion probability for heavier charges. This method fails 
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FIG. 1. The time evolution of (a) mean nucleonic density 
p avg / p a ; (b) size of the heaviest fragment A max ; and (c) the 
multiplicity of clusters with mass A > 5. We display the 
reaction of Au(35 MeV/A)+Au at reduced impact parameter 
b/b max =0.55; b max = 1.142[A T 1/3 + A p 1/z ]. 

to break-up the spectator matter into smaller fragments, 
therefore leading to overestimated charge of heavier frag- 
ments even at 300 fm/c. As a result, the tail of charge 
spectrum shifts towards the higher Z values with MST 
method. This also highlights the discrepancy in MST 
procedure to describe the spectator fragmentation us- 
ing semi-classical transport method. The origin of the 
fragments with SACA method is, however, quite earlier 
determined when nuclear matter is still excited and com- 
pact. 

To derive quantitative information on the phenomenon 
of stopping, one needs to study distribution of frag- 
ments in velocity space. We display in Fig. 3, the 
longitudinal rapidity distributions for light charged par- 
ticles LCPs [2 < A < 4] and intermediate mass frag- 
ments IMFs [5 < A < 65] emitted in reaction of Au(35 
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FIG. 2. The charge dispersion of the nuclear fragments in pe- 
ripheral Au (35 MeV/A)+ Au collisions in the impact param- 
eter interval 0.5 < b/b max < 0.6 (top) and 0.6 < b/b max < 0.7 
(bottom). The model calculations using the MST (dotted 
curve, 300 fm/c) and SACA (solid curve, 100 fm/c) ap- 
proaches are compared with the experimental data (filled cir- 
cles) PI- 

MeV/A)+Au at b/b max = 0.55. Interestingly, we find 
that the MST method predicts the fragment-emission 
from mid-rapidity zone only (as indicated by a peak 
around (y/ybeam)cm = 0). SACA approach, on the other 
hand, indicates significant contribution coming from the 
target-like and projectile-like remnants apart from the 
mid- velocity emission. MST method can't cause cause ef- 
fective break-up of spectator components to generate the 
IMFs close to the target and projectile rapidities even at 
peripheral geometries. This reflects non-equilibrium sit- 
uation in the fragmenting system with SACA approach. 
The emission of the LCPs and IMFs near the projectile 
and target velocities reflects essentially the binary char- 
acter of the collisions, apart from the midrapidity source 
emission. The peaks observed at the target and projec- 
tile rapidities is also indicative of dynamical scenario of 
multi-fragment emission where system has not enough 
time to pass through a state of thermodynamical equi- 
librium. It is worth mentioning that Ca+Ca collisions at 
35 MeV/A were studied recently within antisymmetrized 
molecular dynamics (AMD) model [25j |. Fragment ob- 
servables such as transverse kinetic energy of fragments 
and radial size of the reaction system are observed to 
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FIG. 3. The longitudinal rapidity distribution dNFrag/dy of 
LCPs (top panel) and IMFs (bottom panel) observed in the 
reactions of Au (35 MeV/A)+Au at a reduced impact param- 
eter of b/b m ax = 0.55. Our calculations are performed within 
MST (dotted curve) and SACA (solid curve) approaches. 
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deviate appreciably from the assumption of equilibrium 
ensemble. This study also favored dynamical scenario of 
fragmentation. 

Next, we turn to the cluster distribution obtained us- 
ing MST and SACA approaches in three-dimensional (3- 
D) coordinate space. Figure 4 displays the 3-D snap- 
shots of cluster distribution obtained in a single event 
of Au(35 MeV/A) + Au collision at b/b max = 0.55. In 
the MST method, free nucleons and LCPs [2 < A < 4] 
are abundantly scattered in the whole space (shown as 
red spheres) indicating their isotropic emission from the 
participant zone. Only a small fraction of intermediate 
mass fragment IMFs [5 < A < 65] can be seen coming 
out of the central overlap region. On the other hand, 
a significant enhancement in the production of IMFs is 
clearly visible for spectator zones using the SACA pic- 
ture. The contribution towards the IMFs doesn't seem to 
come from any specific region. In other words, the QMD 
+ SACA calculations suggest that IMFs originate from 



FIG. 4. (Colored online) 3-D snapshots of a single event of 
Au (35 MeV/A) + Au at b/b max = 0.55 using MST (top) 
and SACA (bottom) pictures. Black spheres represent the 
nucleons distribution in IMFs and red spheres represent the 
nucleons distribution in clusters with mass A = 1 — 4. 



the 'extended neck' region as well as from the spectator 
zones. The MST approach fails to break the spectator 
components efficiently, thereby under-estimate the IMF 
yields. These results highlight the importance of clus- 
tering criterion in describing reaction mechanism in low 
energy domain of HI reactions. 

Taking the advantage that SACA approach is able 
to reproduce the experimental charge yields quite accu- 
rately, we further compare mean charge of the heaviest 
fragment (Z max ) obtained using MST and SACA meth- 
ods with data. Figure 5 shows the mean charge of heavi- 
est fragment (Z max ) as a function of 'reduced' impact pa- 
rameter b/b max for Au (35 MeV/A)+Au reactions. The 
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FIG. 5. The charge of the heaviest fragment (Z max ) as 
a function of the reduced impact parameter b/b ma x- Our 
calculations employing the MST (dotted curve) and SACA 
(solid curve) approaches are compared with experimental data 
(filled circles) [la |. 

model predictions using SACA (at 100 fm/c) and MST 
(at 300 fm/c) approaches are displayed along with ex- 
perimental data taken with combined Multics-Miniball 
(MM) array [llj]. One can see an increasing trend of 
r?max w j^jj impact parameter as is expected due to in- 
crease in size of spectator zone in the exit channel. Here 
also, QMD+MST approach fails to break the spectator 
matter effectively, and therefore leads to overpredicted 
Z max even at fm / a The SACA method, on other 

hand, reproduces the experimental trend quite accurately 
at an earlier time. From these finding, we can infer that 
the SACA method is on reliable footing to explore the dy- 
namics of spectator matter fragmentation in low energy 



heavy-ion collisions. Further it allows faster recognition 
of clusters on reaction time scale. 



IV. SUMMARY 

Summarizing, we have performed a comparative study 
of two different clusterization models by simulating 
the peripheral 197 Au+ 197 Au collisions at 35 MeV/A. 
Our calculations within QMD approach coupled with 
SACA clusterization subroutine showed that spectator 
zones contribute significantly towards fragment pro- 
duction in peripheral HI collisions. Contrary to this, 
the MST method could predict fragment-emission from 
mid- velocity source only even at 300 fm/c after the 
initial contact between projectile and target nuclei. 
This questions the validity of employing MST method 
of clusterization to investigate reaction dynamics at 
low incident energies. Our model predictions using 
the SACA method reveal that significant contribution 
towards IMFs emission is located near the initial target 
and projectile rapidities. Our model predictions for 
the charge distribution and mean charge of heaviest 
fragment (Z max ) using the SACA method are in nice 
agreement with experimental data taken with Multics- 
Miniball (MM) array [13, El- These findings highlight 
the importance of clustering criterion in describing 
mechanism behind low-energy spectator fragmentation. 
The SACA formalism can, therefore, be treated as 
general improvement over conventional clusterization 
algorithm to recognize the fragment structure at low 
and intermediate energies. 
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